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ABSTRACT 

The spectrophotometric calibration of surveys is a significant, but often neglected, issue when 
measuring the history of star formation by combining spectroscopic surveys conducted with differ- 
ent instruments. We describe techniques for photometric calibration of optical spectra obtained 
with the MMT's fiber-fed spectrograph, Hectospec. The atmospheric dispersion compensation 
prisms built into the MMT's f/5 wide field corrector effectively eliminate errors due to differ- 
ential refraction, and simplify the calibration procedure. The procedures that we describe here 
are applicable to all 220,000+ spectra obtained to date with Hectospec because the instrument 
response is stable. We estimate the internal error in the Hectospec measurements by comparing 
duplicate measurements of ~1500 galaxies. For a sample of 400 galaxies in the Smithsonian Hec- 
tospec Lensing Survey (SHELS) with a median z=0.10, we compare line and continuum fluxes 
measured by Hectospec through a 1.5" diameter optical fiber with those measured by the Sloan 
Digital Sky Survey (SDSS) through a 3" diameter optical fiber. Agreement of the [Oil] and Ha 
SHELS and SDSS line fluxes, after scaling by the R band flux in the different apertures, suggests 
that the spatial variation in star formation rates over a 1.5 to 3 kpc radial scale is small. The 
median ratio of the Hectospec and SDSS spectra, smoothed over 100 A scales, is remarkably 
constant to ~5% over the range of 3850 to 8000 A. Offsets in the ratio of the median [Oil] and 
Ha fluxes, the equivalent width of US and the continuum index d4000 are a few percent, small 
compared with other sources of scatter. We also explore the impact of atmospheric absorption. 
Observing redwards of 6500 A, it is impossible to remove the effects of atmospheric absorption 
perfectly because the variation of absorption with wavelength is not resolved by a moderate 
dispersion spectrograph. Thus measurements of spectral line fluxes including Ha, and derived 
physical quantities including star formation rates, may have sizable systematic errors where the 
redshifted spectral features land on strong atmospheric absorption troughs. 

Subject headings: techniques: spectroscopic, galaxies: absorption lines, galaxies: emission lines 
Introduction mation rate densities from a variety of surveys 



As spectroscopic and photometric surveys in- 
crease in quality and size, the physical parameters 
which govern galaxy evolution become increas- 
ingly accessible. A range of measures of star for- 



tion history of the universe (e.g. iHoDkins 2004: 


Lv et al. 2007; Tresse et al. 200/1 


Reddv et al. 


2008: Shiova et all 2008: Villar et al. 


2008). 



One of the many issues in reconstructing the 
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star formation history is the relative calibration 
of surveys which use different measures (Hopkins 
2004). After calibration, the scatter in the mea- 
sured star formation rate density at fixed red- 
shift remains large. Even for the same star for- 
mation indicator, inadequacy of the relative cal- 
ibration of surveys with different instruments on 
different telescopes may be one source of scatter 
and systematic offset. Here we use the overlap 
of two large spectroscopic samples from different 
telescopes with different fiber instruments to eval- 
uate these issues for emission line fluxes, US equiv- 
alent width, and the strength of the 4000 A break. 
We demonstrate that with attention to technical 
detail, the median offsets in all of these measures 
are at the few percent level. 

One of our overlapping samples is the SHELS 
survey (IGeller et al.l I2005B carried out with the 



Hectospec (jFabricant et al.l 120051 ) on the MMT; 
the oth er is derived from the Sloan Di gital Sky 
Survey (jAdelman-McCarthv et al1l2008l ). We de- 
scribe the spectrophotometric calibration of Hec- 
tospec and estimate the internal errors in this cal- 
ibration by examining 1467 unique pairs of mea- 
surements of a subset of the galaxies. We estimate 
the external errors in the spectrophotometric cal- 
ibration by comparisons with galaxies observed in 
common with the Sloan Digital S ky Survey (SDSS, 
Adelman-McCarthv et"aD)2008l ). 

We extract the absolute calibration for our 
spectrophotometry from external photometry; we 
use the Deep Lensing Su r vey R -band photome- 



try (jwittmaneiaD Hqq3) ■ This absolute 



calibration corrects for clouds, seeing, telescope 
tracking and guiding, as well as errors in astrome- 
try and alignment, as long as these losses of light 
are independent of wavelength. The spectrophoto- 
metric errors asso ciated with differential refraction 
( Filippenko 19821 ) are largely eli minated because 



the MMT's wide-field corrector ( Fabricant et al 



20041) contains atmospheric dispersion compensa- 
tion (ADC) prisms. The atmospheric dispersion 
across Hectospec's wavelength range (~3700 to 
8500 A), if not corrected by the ADC prisms, 
would exceed the Hectospec's 1.5" fiber diameter 
at a zenith angle of 45°, or an airmass of 1.4. 

Section 2 describes the data and the procedures 
we use to calibrate the Hectospec spectrophotom- 
etry Because Hectospec's throughput is stable 
(Section 2.3), this calibration procedure is appli- 



cable to all Hectospec observations. In Section 3 
we measure the internal errors in the Hectospec 
spectrophotometry by comparing repeated obser- 
vations; in Section 4 we estimate the external er- 
rors in our spectrophotometry by comparing Hec- 
tospec and SDSS spectra. In Section 5 we discuss 
the impact of calibration issues on determination 
of the star formation history. We conclude in Sec- 
tion 6. 

2. Spectrophotometry with Hectospec 

2.1. Introduction 

In this section we describe the procedures we 
use to calibrate the Hectospec spectrophotometry. 
A w eb page, tdc-www. harvard.edu/instruments/hectospec/reduce.htn. 
and lMink et al.l (|2007l ) describe the s tandard CfA 
Hecto spec reduction pipeline (see also lFabricant et aD 
(|2005l )). We further process the spectra emerging 
from this pipeline to convert observed counts to 
flux. There are four steps in this process: (1) we 
correct for the smooth component of atmospheric 
extinction, (2) we remove, as far as possible, the 
sharp H2O and O2 atmospheric absorption lines at 
wavelengths redder than ^6000 A, (3) we correct 
for Hectospec relative throughput as a function 
of wavelength, and (4) we compare a synthetic 
R-band flux derived from the spectra to a corre- 
sponding R band aperture magnitude to deter- 
mine the absolute flux normalization. 

2.2. Data Sources 

To calibrate and test our spectrophotometric 
reductions we use data from three surveys: the 
Smithsonian Hecto spec Lensing Survey (SHELS; 
I Geller et al.l (120051)): the De ep Lens Survey (DLS: 
Iwittman et al.l 120061 . 120021): and the Sloan Digi- 
tal S ky Survey (SDSS: lAdelman-McCarthv et all 
20081 ). 

The SHELS F2 galaxy sample includes spec- 
tra for 9825 galaxies to a limiting total R magni- 
tude of 20.3 in a field of 4 square degrees. The 
R band magnitude s are on the Cousins system, 



Bessell et al.l (|1998[ ). The integral SHELS com- 
pleteness to the limiting magnitude is 97.7%, and 
the differential completeness at the limiting mag- 
nitude is 95.5%. The F2 field covers the RA range 
9.24667 - 9.40389 hr and the DEC range 28.97722° 
- 31.02333°. The spectra were obtained with Hec- 
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tospec on the MMT between April 13, 2004 and 
April 20, 2007. 

From the SHELS dataset we use 2934 SHELS 
spectra (1467 pairs) to calibrate our internal errors 
and 406 spectra for galaxies which also have SPSS 
spectrophotometry (jAdelman-McCarthv et al.l200 
to calibrate our external errors. 

The SHELS galaxies are selected based on DLS 
total R magnitudes; here we use DLS R aperture 
photometry to calibrate our spectra to absolute 
units of flux. We use 1.5" diameter aperture pho- 
tometry to calibrate the flux in our 1.5" fibers, 
and we use 3" aperture magnitudes to transform 
our spectrophotometry to match the spectropho- 
tometry of the SDSS, transformed to their 3" fiber 
apertures. 

The SHELS spectra were obtained with Hec- 
tospec's 270 line mm -1 grating that provides a 
dispersion of 1.2 A pixel -1 and a resolution of 6.2 
A FWHM. Wavelength calibration is carried out 
with a hollow Fe cathode NeAr lamp illuminating 
a screen mounted on the MMT enclosure; typical 
scatter in the fits of the wavelengths of line cen- 
ters to low order polynomials is 0.06 A. Exposures 
range between 45 and 120 minutes. Hectospec's 
300 fibers can be placed within a circular field of 
view 1° in diameter. Typically 10 to 15% of the 
fibers are used for sky background measurement. 
Sky subtraction accuracy is ultimately limited by 
the fiber to fiber throughput variations (2% RMS) 
that are calibrated with dome flats or twilight 
flats. System throughput (Hectospec plus MMT) 
for light entering a 1.5" aperture peaks at 17% at 
-6000 A, falling to half this val ue at 3900 and 



7900 A. Full details are given in iFabricant et al 
(120051) 



2.3. Atmospheric Extinction 

High precision photometry or spectrophotom- 
etry requires frequent measurements of atmo- 
spheric extinction, but the necessary observations 
are costly in telescope time. Instead, we use a 
standard extinction curve obtained at KPNO and 
distributed with IRAF. We do not have a mea- 
sure of the night to night extinction variations, 
and such variation is included when we measure 
internal or external errors in our measurements in 
Sections 3 and 4. The standard extinction curve 
only includes the continuum extinction and not 



the narrow H2O and O2 absorption features in- 
creasingly prominent in the red. The absorption 
from these narrow features can only be removed in 
an average sense from moderate dispersion spec- 
tra; the intensities of spectral lines in the target 
object that fall in the absorption troughs are cor- 
respondingly uncertain. 

Figure [1] shows the narrow absorption features 
between 6500 and 9000 A plotted at 2 A res- 
olution; we also show a spectrum of the night 
sky at the same resolution. The absorption fea- 
tures are taken from a high dispersion atlas ob- 
tained with the Fourier Transform Spectrograph 
at the McMath Solar Observatory available at 
ftp.noao.edu / catalogs /atmospheric_transmission / . 
The upper axis shows the redshift of Ha at the 
wavelengths plotted on the lower axis. Regions 
free of strong night sky emission lines are not nec- 
essarily free of strong absorption features. 

The McMath measurement of atmospheric ab- 
sorption lines is approximately normalized to one 
airmass at an altitude of 2096 meters. At the 
MMT altitude of 2600 meters, the absorption is 
somewhat less. We have empirically renormalized 
the McMath absorption to produce the best av- 
erage correction of the Hectospec spectra for the 
strong absorption features between 6700 and 8000 
A; reducing the McMath absorption by 15% pro- 
vides the best result. 

The standard Hectospec pipeline applies a par- 
tial correction for the atmospheric absorption fea- 
tures, including only the strong features at ^6800 
and ^7600 A. We remove this pipeline correction 
and apply a correction using all of the features in 
the McMath spectrum, binned to a resolution of 
1 A. Because the absorption varies on a wave- 
length scale finer than the resolution of a mod- 
erate dispersion spectrograph, the absorption cor- 
rection is approximate. Although the corrected 
continuum looks smooth, the fluxes from narrow 
emission lines are not measured to correspondingly 
high accuracy. 

We can estimate the potential errors in line 
fluxes for narrow emission lines by calculating the 
difference between the McMath absorption at 0.1 
and 2 A resolution, normalized by the absorption 
at 2 A resolution. For regions with strong absorp- 
tion, for example 6800-6900 AA and 8100-8300 A, 
the mean normalized difference is 7-8%, and the 
RMS normalized difference is 14-17%. 
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2.4. Hectospec Relative Throughput 

We observe spectrophotometric standard stars 
periodically with Hectospec to monitor its through- 
put. In the early operation of the instrument we 
allowed many months between standard star ob- 
servations, although the current practice is to ob- 
serve standard stars at least once per run. Fortu- 
nately, Hectospec's long term stability has proven 
to be excellent as discussed below. No secondary 
spectrophotometric standards were observed si- 
mult aneously with the SHELS galaxy observations 
(e.g. IStoughton et al.ll2002h . 

One complication in interpreting the standard 
star spectra is the presence of a significant amount 
of second order light at wavelengths redwards of 
7000 A when observing blue spectra. We do not 
currently use an order blocking filter with Hec- 
tospec because the loss of spectrophotometric pre- 
cision is outweighed by the desire for maximum 
throughput. 

For maximum spectrophotometric accuracy it 
is necessary to use a standard star with a spec- 
trum as similar as possible to the observed ob- 
jects. Most of the observed spectrophotometric 
standard stars are significantly bluer than the av- 
erage gala xy spectrum in the SHELS survey. One 
standard ijMassev et al.lll988l ). VI Cyg (Cyg OB2 
No. 9), a reddened 051 supergiant, has a very red 
spectrum that reduces second order contamina- 
tion from the grating to a negligible level, and is 
thus well matched to a typical SHELS galaxy. We 
use data obtained from VI Cyg in April 2007 to 
calibrate the Hectospec throughput for the spec- 
trophotometry described here. 

We have spectrophotometric observations of 
BD+284241, a subdwarf O star with a blue spec- 
trum, from widely separated times. These spectra 
allow us to look for changes in instrument through- 
put over this time. Figure [2] shows the ratio of 
spectra obtained in October 2004 and November 
2007, normalized to unity at 6000 A. Differences 
in seeing and extinction complicate this compar- 
ison, but the relative throughput appears to be 
constant to ±10% between 4000 and 8000 A over 
three years. 



2.5. Absolute Spectrophotometric Cali- 
bration 

Prior to the absolute spectrophotometric cal- 
ibration, the spectra processed by the standard 
pipeline are corrected for extinction, converted 
from counts to flux in f\ units (ergs cm -2 sec -1 
A -1 ), and corrected for relative throughput as a 
function of wavelength. The absolute flux nor- 
malization is arbitrary because the observations 
were made in variable seeing and possibly through 
clouds. 

We begin our absolute spectrophometric flux 
calibration by calculating a synthetic R-band flux 
in f\ units. 



jF(X)T(X)dX 
jT(X)d\ 



(1) 



Here, fa is the flux at the R band effective 
wavelength (~6410 A for the Cousins system, 
Bessell et al.1 (|l998h h F(X) is the raw Hectospec 
flux as a function of wavelength, and a is a normal- 
ization constant required to convert the raw Hec- 
tospec flux to a calibrated flux. All fluxes are in i\ 
units. We calculate T(X), the R-ban d response for 



the K PNO Mosaic camera system (jMuller et al 



1998! ), as the product of the K1004 filter response 
and the throughput of the Mosaic CCDs. The final 
step is to determine the correct normalization, a, 
by converting the R-band 1.5" diameter aperture 
magnitudes to an R-band flux (ergs cm -2 sec -1 
A" 1 ). We use the results tabulated in Table A2 of 
Bessell et al.l (|l998l ) (note that the f\ and i v zero 
points are transposed), where 



f R = io(- - 4 *( fl + 21 - 655 )) 



(2) 



2.6. Index Definition 



We use [Oil] A3727, Hi5, and Ha line fluxes and 
equivalent widths, and d4000, a measure of the 
strength of the 4000 A spectral break. For the line 
fluxes and equivalent widths we use indices rather 
than fits to the line profiles. For narrow lines the 
results are freer of systematic fitting errors to the 
flat-topped point spread function arising from the 
resolved fiber image. Table 1 defines the wave- 
length regions used to measure the [Oil] A3727, 
H<5, an d Ha line fluxes and equivalent widths. Fol- 
lowing [BaiogheraD (|l999h . we define d4000, the 
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strength of the 4000 A break, as the ratio of flux 
in the 4000-4100 A band to that in the 3850-3950 
A band. 

3. Hectospec Internal Spectrophotometric 
Errors 

We estimate the internal errors in the Hec- 
tospec spectrophotometry from a sample of 1468 
unique pairs of duplicate SHELS spectra. This 
large set of repeat measurements allows us to esti- 
mate our systematic errors accurately. We include 
only spectra obtained when the ADC prisms were 
operating properly (see Table 2). 

We compare multiple flux measurements of 
[Oil] A3727 and Ha, multiple equivalent width 
measurements of US, and multiple measurements 
of d4000. We propagate the statistical errors for 
each pixel of the spectra for the equivalent width 
and line flux measurements. 

Additional systematic errors arise from several 
sources: short term variations in the wavelength 
dependence and amount of extinction, the wave- 
length dependence of seeing variations, variable 
alignment of the fiber with respect to the observed 
galaxy, errors in our calibration of fiber to fiber 
throughput variations, imperfect sky subtraction, 
and the (small) wavelength dependence of cloud 
extinction. We estimate the magnitude of these 
systematic errors by examining the differences be- 
tween repeated measurements. We add a system- 
atic error in quadrature to the statistical errors to 
calculate a total internal error. We determine the 
total internal error in each measurement, <r t , by 
requiring 68% of the measurement differences to 
lie within \J~2ot ■ The systematic error can then be 
easily calculated. 

We first compare 762 repeated measurements of 
the [Oil] A3727 and Ha emission line fluxes; these 
are plotted in Figure [H We only include measure- 
ments where the EW([OII]) or EW(Ha) >3 A, and 
the emission line is significant at >3cr. We deter- 
mine the systematic error to be 18% of the line 
flux. The total RMS scatter in the measurements 
is 22%, and 6.8% of the measurement differences 
are larger than y/2(2at), to be compared with 4.5% 
for a Gaussian error distribution. 

Figure [4] plots 592 repeated measurements of 
Ha emission line fluxes. We determine the sys- 
tematic error to be 18% of the line flux. The 



total RMS scatter in the measurements is 23%, 
and 9.1% of the measurement differences are larger 
than V2(2a t ). 

Figure [5] plots 1011 repeated measurements of 
H<5 equivalent width where the total measurement 
error (statistical and systematic) is 3 A or less. 
We determine the systematic error to be 2.3 times 
the statistical error. The total RMS scatter in the 
measurements is 1.5 A, and 4.4% of the measure- 
ment differences are larger than ^/2(2<j t )- 

Figure [6] plots all 1468 repeated measurements 
of d4000. We determine the systematic error to 
be 4.5% of the index value. The formal statistical 
error for this measurement is very low, so the total 
error is dominated by the systematic error. The 
total RMS scatter in the measurements is 0.086, 
and 9.8% of the measurement differences are larger 
than \/2(2a t ). 

In all cases, the systematic errors dominate the 
statistical errors. We note that the systematic er- 
rors in the repeated measurements of [Oil] and 
Ha line fluxes are both 18% even though these 
lines are widely separated in observed wavelength. 
The d4000 measurements appear quite robust and 
their errors are considerably smaller than the line 
flux errors. The internal error estimates are sum- 
marized in Table 3. 

4. Comparison with SDSS Spectra 

To evaluate the external errors in our measure- 
ments of emission line fluxes, equivalent widths, 
d4000, and continuum shapes, we compare 406 
galaxies within both the SHELS sample and the 
well-calibrated Sloan Digital Sky Survey (SDSS) 
Data Release 6 (DR6). The set of overlapping 
spect ra includes galaxies from the LRG sur- 
vey (|Eisenstein et al.l l200lh ( 25%) with the re- 
mainder from the m ain sample of SDSS galaxies 
(|Strauss et al.ll2002h . 

We reanalyze the fluxed DR6 SDSS spectra in 
the same fashion as the SHELS spectra to ensure 
a precise comparison with the SHELS results. The 
DR6 and the SHELS flux calibrations are indepen- 
dent of one another. We do not use the published 
SDSS measurements of line fluxes or spectral in- 
dices because our goal is to compare the flux spec- 
tra, not analysis techniques. The median redshift 
of the SHELS/SDSS sample is 0.13 for the full 
sample of 406 galaxies; the median redshift is 0.10 
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for the 154 galaxies with Ha emission. Figure [7] 
is the redshift histogram for the galaxies with Ha 
emission, and for the full sample. 

A major subtlety in comparing the SHELS and 
SDSS spectra is the difference in fiber apertures, 
1.5" and 3.0", respectively. We finesse this prob- 
lem by using the DLS photometry taken with 
seeing between 0.82 and 0.91". The DLS pho- 
tometry provides the 1.5" and 3.0" aperture ri- 
band magnitudes we use for relative calibration 
of the SHELS and SDSS spectra; we scale the 
fluxes through the Hectospec and SDSS apertures 
to these R-band magnitudes. This procedure as- 
sumes that the emission line fluxes scale as the 
R-band light. Physical scatter in this comparison 
certainly results from variation in the distribution 
of emission-line regions relative to the continuum 
emission in the galaxies. 

One potential confusion in comparing the 
SHELS and SDSS spectrophotometry is that the 
SDSS DR6 spectra are flux calibrated to the in- 
tegrated light from standard stars, not the por- 
tion of th e light that arrives within a fibe r aper- 



ture (see I Adelman- McCarthy et all liooih . This 



procedure might be convenient for stars, but is 
less so for extended objects like galaxies. We 
have looked at a selection of SDSS stars in the 
SHELS field to derive a correction factor; the total 
(PSF) magnitudes are 0.34 magnitudes brighter 
than the 3.0" aperture magnitudes. This re- 
sult agrees with the mean differ ence reported in 
Adelman-McCarthv etail (|2008l ). 

The seeing at the time of the SHELS Hectospec 
observations is not a factor in the comparison be- 
tween the SHELS and SDSS spectra because the 
SHELS spectrophotometry is normalized to the 
DLS photometry. In fact, the normalization of the 
flux entering the Hectospec 1.5" aperture drops 
out of the SHELS-SDSS comparison because we 
use the R-band photometry to convert the SHELS 
flux to a 3.0" aperture. 

Figure M plots the median ratio of 273 SHELS 
and SDSS galaxy spectra between 3850 and 8200 
A, boxcar smoothed over 100 A intervals. All pairs 
of spectra with at least 25 counts per pixel in the 
raw SHELS spectra over the entire spectral range 
are included in the median. The ratios of the spec- 
tra are normalized to unity between 5000 and 5500 
A. The interquartile range is also plotted. On av- 
erage, the SHELS fluxing procedure preserves the 



spectral shape to high accuracy, on the order of 
±5%. 

The comparison of [Oil] A3727 and Ha line 
fluxes measured from the SHELS and SDSS spec- 
tra shows that the SHELS flux calibration proce- 
dures work remarkably well. Figure [9] plots SDSS 
[Oil] lines fluxes against SHELS line fluxes for 128 
galaxies with [OH] equivalent width >3 A and 
are significant at >3 a confidence. The ratio of 
the SHELS median flux to the SDSS median flux 
is 0.96. We calculate a systematic error for the 
SDSS line fluxes and equivalent widths in same 
fashion as we did above for the SHELS systematic 
errors, but this time using the total errors deter- 
mined above for the SHELS data. The assignment 
of this additional systematic to the SDSS data is 
not above question, but we believe that system- 
atic scatter affecting the SHELS data should have 
shown up in our internal comparisons. The statis- 
tical errors for the SDSS measurements are prop- 
agated from the errors associated with each pixel. 

We determine a systematic error of 12% for the 
SDSS [Oil] flux measurements. The total RMS 
scatter in the [Oil] flux measurements is 22% (as- 
suming an equal division of errors between SDSS 
and SHELS; the RMS scatter in the differences be- 
tween the measurements is larger). The frac- 
tion of >2<t outliers is 4.7%. 

Figure [TU] is the corresponding plot for the Ha 
line fluxes of 154 galaxies where the line has an 
equivalent width >3 A and is significant at >3 
a. Here, the ratio of median SHELS to SDSS line 
flux is 1.03, and an additional systematic error of 
15% is determined for the SDSS line fluxes. The 
total RMS scatter for measurements is 22%, and 
the fraction of >2a outliers is 9.1%. 

Figure [Tl] compares the H5 equivalent widths 
measured from the SHELS and SDSS spectra. We 
plot data for 352 galaxies with total measure- 
ment errors (statistical and systematic) of 3 A or 
less. The ratio of the SHELS to the SDSS median 
equivalent width is 0.97. We calculate that the 
SDSS total error is 1.2 times the statistical error. 
The total RMS scatter per measurement is 1.0 A, 
and the fraction of >2a outliers is 5.1%. 

Figure [T2] compares the d4000 measured from 
SDSS and SHELS spectra for 358 galaxies. The 
ratio of the median SHELS to the median SDSS 
d4000 measurement is 1.0. We determine that the 
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SDSS total error in d4000 is 0.03 times the index 
value. The total RMS scatter per measurement is 
0.064, and the fraction of >2cr outliers is 5.6%. 

Table 4 summarizes the comparison of the 
SHELS and SDSS measurements. Our compar- 
isons of emission line fluxes and continuum fea- 
tures in spectra of the same objects taken from the 
SDSS DR6 and SHELS show that spectrophoto- 
metric measurements from these two fiber surveys 
agree remarkably well; there are no significant 
offsets. The calibration procedures for the two 
instruments differ but yield consistent results. 

The errors in emission-line fluxes substantially 
exceed errors in the continuum measure, d4000. 
This difference results in part from both the dis- 
creteness of the emission-line regions and the vari- 
ation in their distribution relative to the contin- 
uum light we use to normalize the SHELS spectra. 
The very small errors in d4000 and the absence of 
any detectable variation relative to the larger aper- 
ture SDSS photometry suggest that the underly- 
ing stellar population has little variation over the 
radial scales, 1.4 and 2.8 kpc (at z=0.1), sampled 
by the apertures. The small d4000 measurement 
errors and its comparative immunity to aperture 
effects encourage its us e as a robust measure of the 
star formation history (jKauffmann et al. 2003). 

5. Astrophysical Implications 
5.1. Introduction 

Understanding star formation in galaxies and 
galaxy evolution re quires spectrophotometry (e.g. 
Kennicutt 1992aHbl:lGavazzi et al. 2002; Shiova et al.l 



lution issues can be profitably explored from red- 
shift zero to 0.7 or more; for quasar spectroscopy 
the redshift range is much larger. As examples of 
the possible impact of combining such surveys, we 
comment on the impact of aperture effects and on 
some issues in the determination of star formation 
rate densities. 

5.2. Aperture Effects 

Aperture effects are an important limitation in 
comparing samples of galaxies over a broad red- 
shift range. iKewlev et al. ( 2005 ) discuss the com- 
plex range of issues which result from a spec- 
troscopic aperture of fixed angular size rather 
than a fixed physical size. By comparing nuclear 
spectroscopy with integrated spectr oscopy of the 
Near by Field Galaxy Sample (NFGS: l"jansen et al 



2000), they demonstrate that apertures enclosing 



20% or more of the galaxy light provide an unbi- 
ased view of the star formation properties. They 
emphasize that the physical aperture which in- 
cludes 20% of the light is a function not only of 
galaxy type but also of galaxy luminosity. Thus 
at greater redshift where surveys are limited to 
the intrinsically most luminous objects, a larger 
physical aperture is necessary. 

( 2005h 



Figures 6 and 7 of IKewlev et al 



are 

useful for judging the adequacy of the Hectospec 
and SDSS apertures as a function of redshift and 
galaxy properties. They show that for the mean 
NFGS galaxy, the SDSS fibers include more than 
20% of the light at z>0.05; the Hectospec fibers 
cover this fraction at z>0.1. For galaxi es brighter 
than t he characteristic luminosity in the lSchechter 



2002|lLamareille et all2005|]2006tlMoustakas fc Kennicuttl(|l976l ) parameterization, L*, the SDSS fibers in 
20061 iMoustakas et alj 120061 : IWeiner et all l2007t 
Cooper et al.ll2008f) . Star formation rates and red- 



dening estimates require line flux measurements. 
Exploring these astrophysical issues across a large 
redshift range generally requires a combination of 
surveys from different telescopes. These compar- 
isons may be compromised without a demonstra- 
tion that the surveys are properly calibrated with 
respect to one another. 

With attention to a range of technical details, 
fiber spectroscopy provides robust line fluxes. In 
the case of Hectospec, spectrophotometry in the 
fiber aperture is possible because the instrument 
includes ADC prisms. By combining SDSS spec- 
troscopy with Hectospec data, many galaxy evo- 



clude at least 20% of the light only for z>0.1 and 
the Hectospec correspondingly at z>0.2. 

The good agreement of the SHELS and SDSS 
[OH] A3727 and Ha line fluxes indicates that star 
formation rates within apertures of 1.5" and 3" 
scale with the R-band light. The median line 
fluxes agree within 3-4% and, for high signal-to- 
noise spectra, the scatter is dominated by system- 
atics which imply a ~18% typical error. 

At z=0.1, the Hectospec and SDSS apertures 
correspond to radii of 1.4 and 2.8 kpc, respectively. 
The median 25% light radius for the sample galax- 
ies is 2.1 kpc. Thus the SDSS aperture generally 
includes 20% of the light, but the Hectospec aper- 
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ture includes less than the iKewlev et al. (|2005l) 
fiducial fraction. 

Although the median luminosity of the sample 
galaxies is R = -21.5, approximately L* (Brown et 
al. 2001), and the apertures are small, aperture 
effects appear negligible in this dataset. Figure [T3l 
shows the ratio of Ha fluxes for the SDSS and 
SHELS as a function of the absolute magnitude of 
the objects. Even though the SHELS fibers gen- 
erally include less than 20% of the light, there is 
no apparent dependence of the flux ratio on the 
absolute magnitude. We have also checked that 
the flux ratio is independent of redshift, apparent 
magnitude, and 25% light radius. 

5.3. Star Formation Rates, Reddening, 
and Star Formation History 

Star formation rate indicators are a cornerstone 
of studies of galaxy evolution. We focus on issues 
underlying the use of optical spectroscopic indica- 
tors Ha and [Oil]. 

A variety of star formation studies use equiva- 
lent widths of these lines as star formation indica - 
tors (IPoggianti et al.lll999HMouhcine et al.ll2005l ). 
A drawback of this procedure is the lack of red- 
dening and metallicity me a surem ents that require 



line fluxes. IKewlev et al. (2004) discuss some of 



the pitfalls of assuming a reddening independent 
of galaxy properties and redshift. An additional 
issue not widely discussed in the literature is the 
systematic impact of the atmospheric absorption 
features on Ha (see Figure [lj; these effects are a 
complex function of redshift. 

Measures of reddening provided by spectroscopy 
are also important for surveys covering a large 
redshift range. In general, surveys include intrin- 
sically more luminous objects at greater redshift; 
these objects tend to be dustier and more metal 
rich (|Calzettill2001r ). Thus the common assump- 
tions of fixed reddening and metallicity can bias 
our view of the star formation history. Our anal- 
ysis indicates that for high signal-to-noise spectra 
where the systematic error dominates, the la error 
in Ha/H/5 is -25%. 

Strong atmospheric absorption is an additional 
serious issue at redshifts near 0.05, 0.1, 0.16, and 
0.25. The fact that the SHELS and SDSS Ha line 
fluxes agree well does not demonstrate that the 
narrow atmospheric absorption features are well 



corrected. Both surveys lack the spectral resolu- 
tion to remove the effect of these features with 
standard techniques. 

Figure [1] shows that the corresponding wave- 
length ranges are relatively free of strong night sky 
emission, but they fall in bands of significant at- 
mospheric absorption. A number of imaging sur- 
veys geared toward measuring the Ha luminos- 
ity function as a route toward a star formation 
rate density may be biased by the im pact of these 
absorp tion features. For example, IShiova et al 



(|2008f ) make a state-of-the-art measurement of the 
Ha luminosity function for the COSMOS survey. 
They use a narrow band filter centered at 8150 
A with a 120 A width. Their survey is thus cen- 
tered at z = 0.24. Figure Q] shows that at many 
redshifts within the bandpass the atmospheric ab- 
sorption is 20 - 30%. This systematic bias toward 
underestimating the volume averaged star forma- 
tion rate is unavoidable in an imaging survey that 
overlaps the atmospheric absorption features. Un- 
less the redshifts of individual objects are known 
from spectroscopic measurement, it is impossible 
to make a correction for absorption a posteriori 
because the impact of these features is a com- 
plex function of the redshift. Without redshifts, 
it is not possible to estimate overall impact of the 
bias. These absorption features also impact [Oil] 
and Lyq imaging surveys in the same bandp ass 
(jMuravama et alJl2007t iTakahashi et al.ll2007l ). 

We conclude that an unbiased measurement of 
the volume average star formation rate from an 
imaging survey requires knowledge of the redshifts 
of the objects. Spectroscopy enables an approxi- 
mate correction for the systematic underestimate 
of line fluxes. It also provides an independent esti- 
mate of the completeness of catalogs derived from 
the narrow-band images. 

6. Conclusions 

The most difficult spectrophotometric issues 
hampering astrophysical measurements such as 
the evolution of the star formation rate density 
affect all instruments, not just those with opti- 
cal fibers. These issues include: (1) calibration 
across different surveys and several instruments, 
(2) aperture effects on the derived star formation 
rate, and (3) the highly structured atmospheric 
absorption bands redwards of 6000 A. 
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Because of Hectospec's stability and the MMT's 
f/5 corrector with ADC prisms that eliminate er- 
rors from differential refraction, a relatively simple 
procedure for flux calibration of Hectospec spec- 
tra, using infrequent observations of spectropho- 
tometric flux standards for calibration of relative 
throughput and reference to R-band photome- 
try for absolute throughput, works surprisingly 
well. Spectrophotometric calibration of instru- 
ments without ADC prisms requires observation 
of secondary spectrophotometric standards in each 
field. The success of the procedure we describe, 
when tested against well-calibrated SDSS spectra, 
demonstrates that short term extinction varia- 
tions are not a serious issue. Over its 4.5 year 
operational lifetime, the Hectospec has acquired 
over 220,000 spectra that can be calibrated using 
our procedure. 

We recover the SDSS continuum shapes and 
line fluxes to ~ ±10% for a subset of the SHELS 
galaxies also observed by the SDSS, allowing only 
a rescaling between the different aperture sizes 
based on R-band aperture magnitudes. The me- 
dian emission line fluxes, H<5 equivalent widths, 
and the amplitude of the continuum break, d4000, 
for the SDSS and SHELS spectra agree to within 
a few percent. For high signal-to- noise SHELS 
spectra the typical errors in emission line fluxes 
are 18%. We plan to use these measures to fur- 
ther explore spectroscopic investigation of the star 
formation history over the range of redshifts cov- 
ered by the combin ed SDSS and SHELS surveys 
(jKewlev et al.ll2009h . 
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Fig. 1. — Spectrum of sky absorption features and 
of the night sky plotted at 2 A resolution. The 
night sky spectrum is in units of relative f„. 
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Fig. 3. — Comparison of 762 repeated observa- 
tions of [Oil] A3727 line fluxes. We include only 
measurements where the EW [Oil] >3 A and >3 a 
significance. Representative error bars are shown. 
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Fig. 2. — Ratio of two spectrophotometric stan- 
dard star observations (BD+284241) in October 
2004 and November 2007. 
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Fig. 4. — Comparison of 592 repeated observations 
of Ha line fluxes. We include measurements with 
EW Ha >3 A and >3 a significance. 
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H(5 EW (A) 

Fig. 5. — Comparison of 1011 repeated observa- 
tions of US equivalent width. Measurements with 
a total error <3 A arc included. 




Fig. 7. — Histogram of redshifts in the full sam- 
ple of 406 galaxies observed by both SHELS and 
SDSS (light line). Histogram of redshifts for the 
154 galaxies used to compare SHELS and SDSS 
Ha line fluxes (heavy line). The median redshift 
for the full sample is 0.13; the median redshift for 
the Ha line flux sample is 0.10. 
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Fig. 6. — Comparison of 1468 repeated observa- 
tions of d4000, the ratio of flux in the in the 4000- 
4100 A band to that in the 3850-3950 A band. 
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Fig. 8.— Median ratio of SHELS to SDSS spectra 
for 273 galaxies, and the interquartile range. Only 
those SHELS spectra with least 25 ADU per pixel 
between 3850 to 8000 A arc included. The ratio is 
normalized to one between 5000 and 5500 A. 
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[Oil] Flux SHELS (ergs cm" 2 s" 1 A -1 ) 

Fig. 9. — Comparison of [Oil] A3727 line fluxes for 
128 galaxies observed by SHELS and SDSS. We in- 
clude those measurements where the EW [Oil] >3 
A and >3 a significance. The ratio of the median 
SDSS and SHELS measurements is 0.96. 




Fig. 10. — Comparison of Ha line fluxes for 154 
galaxies observed by SHELS and SDSS. We in- 
clude those measurements where the EW [Oil] >3 
A and >3 a significance. The ratio of the median 
SDSS and SHELS measurements is 1.03. 
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Fig. 11. — Comparison of US equivalent widths for 
352 galaxies observed by SHELS and SDSS. Those 
galaxies with equivalent width measurement errors 
<3 A are included. The ratio of the median SDSS 
and SHELS measurements is 0.97. 
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Fig. 12. — Comparison of d4000 for 358 galaxies 
observed by SHELS and SDSS. The ratio of the 
median SDSS and SHELS measurements is 1.00. 
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Fig. 13.— Ratio of SDSS to SHELS Ra fluxes as a 
function of R absolute magnitude. No systematic 
trend is observed, even though the SHELS fibers 
generally intercept less than 20% of the light from 
these galaxies. 



Table 1 

Definition of Spectral Indices (All given in A) 



Index 


Blue Cont. Center (Width) 


Feature Center (Width) 


Red Cont. Center (Width) 


[Oil] 


3685.0 (25.0) 


3727.3 (8.0) 


3770.0 (25.0) 


R6 A 


4060.7 (19.1) 


4102.9 (19.4) 


4144.8 (16.3) 


Ha 


6510.0 (25.0) 


6562.8 (8.0) 


6622.0 (25.0) 



Table 2 

ADC Prism Malfunction Intervals 



Year 


Interval 


Nature of Error 


Effect on Data 


2004 


before 10/30 


ADC prisms rotating in wrong sense 


serious 


2005 


10/25 


ADC prisms off, no correction 


small for low airmass 


2006 


11/16 to 11/22 


ADC prisms off, unknown orientation 


unknown 


2007 


2/21 to 3/16 


ADC prisms off, unknown orientation 


unknown 



Table 3 

SHELS Spectrophotometry and Spectral Index Internal Error Estimates 



Feature 


Number of Pairs 


SHELS Systematic Error 


Total RMS Scatter 


Fraction >2<r 


[Oil] line flux 


762 


18% of flux 


22% 


6.8% 


Ha line flux 


592 


18% of flux 


23% 


9.1% 


H<5 EW 


1011 


total error is 2.3 


1.5 A 


4.4% 






times statistical error 






d4000 


1468 


total error is 0.045 


0.086 


9.8% 






times value 







Table 4 

Comparison of SHELS and SDSS Spectrophotometry and Spectral Indices 



Feature 


Number of Pairs 


Median Ratio 


SDSS Systematic Error 


Total RMS Scatter 


Fraction >2cr 


[OH] line flux 


128 


0.96 


12% of flux 


22% 


4.7% 


Ha line flux 


154 


1.03 


15% of flux 


22% 


9.1% 


H<5 EW 


352 


0.97 


total error is 1.2 


1.0 A 


5.1% 








times statistical error 






d4000 


358 


1.00 


total error is 0.03 


0.064 


5.6% 








times value 
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